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THE  NAVY  TESSELLATED  SPHEROID  MAP  PROJECTION  SYSTEM: 
A  COMPREHENSIVE  DEFINITION 


BACKGROUND 

With  the  introduction  of  digital  moving  map  systems  aboard  tactical  Navy  and  Marine  Corps 
aircraft,  a  pilot’s  use  of  paper  map  products  is  quickly  changing  to  a  reliance  on  digital  maps.  A 
digital  moving  map  system  provides  a  pilot  with  precise,  hands-off  navigation  and  route  selection. 
This  system  was  originally  designed  for  night-attack  air  missions,  since  it  provides  valuable  positional 
information  at  night  when  visual  flight  navigation  is  often  impossible.  A  digital  map  computer 
stores  all  of  the  charts  that  are  required  by  the  pilot  for  a  particular  mission.  Additionally,  the 
system  allows  the  pilot  to  display  new  threat  and  intelligence  information,  or  other  data,  as  overlays 
onto  the  base  chart. 

The  Map  Data  Formatting  Facility  (MDFF)  at  the  Naval  Research  Laboratory  (NRL)  is  developing 
the  Compressed  Aeronautical  Chart  (CAC)  database  in  direct  support  of  aircraft  digital  moving  map 
systems.  CAC  is  derived  from  Defense  Mapping  Agency  (DMA)  ARC  (Equal  Arc-second  Raster 
Chart)  Digitized  Raster  Graphics  (ADRG),  which  consists  of  full-color  charts  that  have  been  scanned 
into  100-pm  pixels.  Each  pixel  consists  of  a  red,  green,  and  blue  (RGB)  component.  The  MDFF 
“downsamples”  ADRG  data  from  100-pm  pixels  (which  provide  an  image  resolution  of  254  ppi 
[pixels  per  inch])  to  approximately  200-pm  pixels  (exactly  128  ppi),  which  results  in  a  4:1  reduction 
in  data  storage.  The  lower  resolution  data  are  then  color-compressed  by  a  factor  of  3:1  and,  finally, 
spatially  compressed  by  a  factor  of  4:1  (Lohrenz  et  al.  1991).  The  final  CAC  product  represents  a 
data  storage  compression  of  48:1  over  the  original  ADRG.  Readers  are  referred  to  Lohrenz  and 
Ryan  (1990)  for  a  detailed  CAC  product  specification. 

CAC  is  stored  as  a  seamless,  global  database  such  that  the  transition  from  one  source  chart  to 
another  is  transparent  to  the  pilot.  The  data  are  stored  in  discrete  51  mm^  (2  in.^)  sections  of  chart, 
called  tessellations  or  segments,  which  cover  the  entire  globe;  the  model  is  known  as  a  tessellated 
sphere  (Fig.  1).  Approximately  four  chart  segments  may  be  projected  onto  the  digital  moving  map 
display  at  one  time. 

The  purpose  of  this  report  is  to  completely  describe  the  CAC  Tessellated  Spheroid  (TS)*  map 
projection  system  so  it  can  be  applied  to  other  databases  of  interest  to  mission  planners.  The  TS 
model  used  for  the  CAC  database  was  originally  referred  to  as  the  Model  IV  Tessellated  Sphere, 
or  TS-4  (Ronish  1986).  It  has  since  been  shortened  to  TS.  The  TS  model  was  originally  conceived 
by  Honeywell,  Inc.  (Ronish  1986;  Honeywell  1990)  for  the  CAC  database. 


*Note  that,  in  this  report,  ”7$”  refers  to  the  specific  TesseUated  Spheroid  model  that  is  used  for  the  CAC  daubase. 
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INTRODUCTION  TO  TESSELLATED  SPHEROID 

Tessellated  spheroids  provide  a  rectangular  coordinate  and  projection  system  at  any  scale  for 
the  entire  earth  ellipsoid.  During  the  conversion  from  ADRG  to  CAC,  data  are  transformed  from 
one  tessellated  spheroid  model  (the  DMA  ARC  projection  system)  into  TS.  Although  there  are 
several  characteristic  differences  between  the  ARC  and  TS  models,  as  presented  in  Table  1,  the 
models  are  functionally  similar. 

a.  Both  are  rectangular  coordinate  systems  for  the  entire  earth  ellipsoid. 

b.  For  both,  the  earth  is  divided  into  latitudinal  bands  with  a  slight  overlap  between  adjacent 
bands. 

c.  Bands  in  both  systems  are  subdivided  into  segments  that  are  bounded  by  fixed  latitude  and 
longitude  lines. 

d.  Segments  in  both  systems  are  divided  into  pixels.  Each  pixel  is  associated  with  a  fixed 
latitude  and  longitude  point  and  an  RGB  value. 

e.  All  latitude  and  longitude  coordinates  in  both  ADRG  and  CAC  have  been  adjusted  to  con¬ 
form  to  the  World  Geodetic  System  of  1984  (WGS  84)  if  the  source  charts  do  not  already  conform 
to  this  standard. 


Table  1  —  Major  Differences  Between  ARC  and  TS 


Categories 

ARC  Values 
(DMA  1989) 

TS  Values 

#  Latitude  Bands 

Pixel  Spacing 

Segment  Size 

#  Pixels/Segment 

Polar  Zone  Projection 

18 

100  ^m 

12.7  mm  (O.S  in.)  square 
127  X 127 

Polar  Azimuthal 
Equidistant 

5 

200  pm 

SI  mm  (2  in.)  square 
256  X  256 

Rotated 

EQZone 
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TS  ZONE  BOUNDARIES  AND 
OVERLAP  REGIONS 

The  TS  projection  system  divides  the  surface  of 
the  earth  into  Hve  latitudinal  bands,  called  zones, 
which  are  numbered  from  0  to  4.  Two  of  these 
zones — North  Polar  (NP)  and  South  Polar  (SP) — 
cover  the  poles,  and  each  of  the  three  nonpolar 
zones — North  Temperate  (NT),  Equatorial  (EQ),  and 
South  Temperate  (ST) — completely  encircles  the 
earth  between  two  fixed  latitude  limits  (see  Fig.  2 
and  Table  2). 

A  region  of  overlap  is  defined  between  each 
pair  of  neighboring  TS  zones,  as  indicated  in 
Fig.  2  —  TS  zones  with  latitudinal  boundaries  and  Figs.  1  and  2.  This  overlap  consists  of  an  extension 

overlap  definitions  of  one  zone  into  an  adjacent  zone  by  two  rows  of 

segments,  where  each  segment  consists  of  2S6  x  256 
pixels.  Thus,  each  region  of  overlap  consists  of  a  total  of  1024  pixels  (512  pixels  overlap  northward 
from  the  southern  zone,  and  512  pixels  overlap  southward  from  the  northern  zone).  For  example, 
two  additional  rows  of  NT  segments  are  added  to  the  southern  boundary  of  the  NT  zone  (at 
31.385°  N)  and  extend  into  the  EQ  zone.  Similarly,  two  additional  rows  of  EQ  segments  are  added 
to  the  northern  boundary  of  the  EQ  zone  (at  31.385°  N)  and  extend  into  the  NT  zone.  Therefore, 
the  CAC  database  stores  each  section  of  overlap  data  twice:  once  as  part  of  one  zone  and  again  as 
part  of  the  adjacent  zone.  The  amount  of  geographic  overlap  (in  degrees  of  latitude)  between  zones 
is  dependent  upon  the  scale  of  the  data  (see  Table  3). 


TS  COORDINATES 


Coordinates  in  the  nonpolar  zones  are  directly  proportional  to  the  WGS  84  coordinates  0  and 
X  (latitude  and  longitude).  Specifically,  the  latitudin^  and  longitudinal  spacing  of  pixels  is  constant 
throughout  any  given  nonpolar  zone.  Coordinates  in  the  polar  zones  have  been  derived  by  a  rotation 
of  the  equatorial  zone’s  coordinate  grid  over  the  polar  regions,  as  defined  in  a  later  section  of  this 
report.  WGS  84  coordinates  6  and  X  in  CAC  are  signed  values  such  that  -90°  <  6  ^  •<■90°,  and 
-180°  ^  X^  •fl80°.  The  X  values  outside  this  range  should 


be  adjusted  by  ±360°  as  appropriate. 


Table  3  —  Extents  of  TS 


Zone  Overlap 
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TS  SEGMENTS 
Segment  Organization 

Each  TS  zone  is  divided  into  spherical  rectangular  segments.  In  the  equatorial  and  temperate 
zones,  these  segments  are  arranged  by  rows  and  columns,  the  boundaries  of  which  are  lines  of 
constant  latitude  and  longitude.  Figure  3  presents  a  section  of  nonpolar  CAC  data  overlaid  with  a 
grid  of  segments,  and  Fig.  4  presents  a  section  of  polar  CAC  data  overlaid  with  segments.  The  polar 
zone  segment  boundaries  are  defined  by  a  rotation  of  the  equatorial  segment  grid  onto  the  polar  regions. 
This  rotated  equatorial  grid  was  extended  to  cover  the  entire  polar  zone,  which  is  larger  than  the 
equatorial  zone  by  approximately  13.85°  of  latitude.  Figure  5  depicts  the  orientation  of  both  grids 
(a  polar  coordinate  grid  overlaid  by  a  rectilinear,  rotated  equatorial,  segment  grid)  for  the  entire  NP 
region. 

Segment  Numbering 

Nonpolar  segments  are  arranged  by  latitudinally  based  rows  and  longitudinally  based  columns.  Polar 
zone  segments  are  arranged  by  latitudinally  based  rows  and  longitudinally  based  columns  within  the 
rotated  equatorial  grid  system.  In  both  cases,  row  0  is  located  with  its  southern  boundary  on 
the  equator,  as  shown  in  Fig.  6a.  Positively  numbered  rows  extend  northward  from  row  0  in  the 
northern  hemisphere,  and  negatively  numbered  rows  extend  southward  from  row  -1  in  the  southern 
hemisphere.  The  first  column  of  segments  in  each  nonpolar  zone  (column  0)  is  located  with  its 
western  boundary  on  the  0“  meridian  (Fig.  6b).  Positive  columns  extend  eastward  from  column  0 
and  stop  at  the  180“  meridian  (Fig.  6c).  Negative  columns  extend  westward  from  column  -1  and 
stop  at  the  180“  meridian.  Table  4  lists  the  TS  scale  models  and  gives  for  each  the  number  of 
segment  rows  and  columns  in  the  equatorial  and  temperate  zones,  excluding  overlap.  To  include 
overlap,  the  number  of  rows  would  be  increased  by  four  (two  rows  for  overlap  on  the  northern  edge 
of  the  zone,  plus  two  rows  for  overlap  on  the  sou^em  edge).  Calculations  for  the  numbers  of  rows 
and  colunms,  as  given  in  Table  4,  are  provided  later  in  this  report. 

Segment  Files 

After  color  compression,  each  TS  segment  is  stored  (in  the  CAC  database)  as  a  single,  compressed, 
raster  image  file.  The  compressed  segment  file-naming  convention  is  xxxxxxxx.xxz,  where  xxxxxxxx.xx 
is  the  same  key  name  that  was  previously  described  (xxxxxxxxxx  in  the  uncompressed  fllename), 
and  z  defines  the  TS  zone  in  which  the  map  segment  resides  (z  =  a  value  between  0  and  4,  as  listed 
in  Table  2).  The  key  name  (xxxxxxxxxx  in  the  uncompressed  file,  or  xxxxxxxx.xx  in  the  compiessed 
file)  is  derived  as  follows: 

KEY  »  [((2  X  MaxRow)  +  1)  x  (COL  +  MaxRow))  +  ROW  +  MaxRow  +  1 

where  MaxRow  a  9000  and  ROW,  COL  a  row  and  column  numbers  of  TS  segment.  Conversely, 
the  segment’s  row  and  column  numbers  can  be  derived  from  the  segment  filename  key: 

ROW  =  mod  ((KEY  -  162018001),  18001) 

COL  a  ((KEY- 162018001) -ROW)/18001 

where  mod  (short  for  modulus)  is  a  software  function  that  returns  the  quantity  that  remains  after 
the  first  argument  (KEY  -162018001)  is  divided  by  the  second  (18001). 
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Fig.  3  —  Nonpolar  CAC  data  overlaid  with  segment  grid 
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Fig.  4  —  Polai  CAC  data  overlaid  with  segment  grid 
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North  polar  zone  coordinate  system  overlaid  with 
rotated  equatorial  segment  grid 
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Fig.  6  —  TS  segment  row  and  column  numbering  scheme;  (a)  origin  for  row  numbering,  (b)  column 
numbering  at  0"  longitude,  and  (c)  column  numbering  at  ISO”  longitude  (1:2M  scale.  EQ  zone) 


Table  4  —  Numbers  of  Segment  Rows  and  Columns  in  Each  I'S 

Scale  Model 


Equatorial 

Temperate 

Scale 

#  Cols  of 
Segments 

#  Rows  of 
Segments 

#  Cols  of 
Segments 

#  Rows  of 
Segments 

1:2,000.000 

380 

68 

304 

22 

1:1.000,000 

760 

136 

608 

44 

1:500,000 

1520 

272 

1216 

88 

1:250,000 

3040 

544 

2432 

176 

1:100,000 

7600 

1360 

6080 

440 

1:50,000 

15200 

2720 

12160 

880 
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Determining  Segment  Size 

Six  scales  of  chart  are  included  in  the  CAC  database,  ranging  from  1  ;50k  (k  =  thousand)  to 
1:2M  (M  =  million)  (Table  5).  The  scale  of  1:2M  is  used  as  the  base  scale  in  calculating  TS 
segment  sizes;  each  segment  in  the  1:2M  scale  ntodel  contains  an  integral  number  of  segments  of 
any  other  scale  model.  Data  at  any  other  scale  that  will  utilize  the  TS  system  must  be  an  integral 
divisor  of  this  base  scale.  Table  6  lists  the  TS  scale  models  that  are  used  in  the  CAC  database  and 
gives  their  segment  sizes  for  equatorial  and  temperate  zones.  The  following  discussion  explains 
how  these  TS  segment  sizes  are  derived.  For  inexact  values,  excess  significant  digits  are  retained 
to  maximize  accuracy  in  the  calculations.  Unless  otherwise  noted,  parameters  were  first  defined  by 
Ronish  (1986),  and  all  conversions  between  English  and  metric  units  of  measurement  are  taken 
from  Beyer  (1979). 

Four  TS  segments  can  be  displayed  simultaneously  on  the  cockpit  screen,  which  is  127  mm^ 
(5  in.^)  square.  A  “human  factors  number"  of  magnitude  1.316,  which  is  incorporated  into  the  area 
coverage  calculations,  enlarges  the  data  and  performs  an  adjustment  for  the  fact  that  the  cockpit 
screen  can  display  only  484  x  484  pixels  of  the  four  segments’  512  x  512  pixel  area  at  one  time. 
The  length  of  one  edge  of  a  512  x  512  pixel  area  is  (127  mm/1.316)  x  (512/484)  =  102.0874677  mm^ 
(4.02  in.2). 


Table  5  —  TS  Scale  Models  and  Source  Charts 
Used  in  CAC 


Scale 

DMA  Source  Chart 

1:2,000,000 

1:1,000,000 

1:500,000 

1:250,000 

1:100,000 

1:50,000 

JNC  (Jet  Navigation  Chart) 

ONC  (Operational  Navigation  Charts 
TPC  (Tactical  Pilotage  Chart) 

JOG  (Joint  Operational  Graphic) 

TLM  (Topographic  Line  Map) 

TLM  (Topographic  Line  Map) 

Table  6  —  Segment  Sizes  for  Each  TS  Scale  Model 


Scale 

Equatorial  Zone 

Temperate  Zones 

Height 

(“lat) 

Width 

("Ion) 

Height 

("lat) 

Width 

(“Ion) 

1:2,000,000 

0.92307692 

0.94736842 

0.92307692 

1.18421053 

1:1,000,000 

0.46153846 

0.47368421 

0.46153846 

0.59210526 

1:500,000 

0.23076923 

0.23684211 

0.23076923 

0.29605263 

1:250,000 

0.11538462 

0.11842105 

0.11538462 

0.14802632 

1:100,000 

0.04615385 

0.04736842 

0.04615385 

0.05921053 

1:50,000 

0.02307692 

0.02368421 

0.02307692 

0.02960526 
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Therefore,  one  segment  (which  consists  of  256  x  256  pixels)  covers  a  102.0874677  mm^  = 
51.04373385  ram^  (2.01  in.^)  portion  of  a  paper  chart.  The  TS  system  identifies  the  location  of  each 
segment  by  a  row  and  column  number,  which  corresponds  to  latitudinal  (LAT)  and  longitu¬ 
dinal  (LON)  location,  respectively.  The  LAT  and  LON  extent  of  each  segment  is  based  on  an 
attempt  to  approximate  a  51.04373385  mm^  area.  The  radius  and  circumference  of  the  earth  (assuming 
a  spherical  earth),  represented  at  a  scale  of  1;2M,  are  required  in  the  row  and  column  calculations. 
Ronish  (1986)  assumed  an  earth  radius  of  125.3264963  in.  at  the  1:2M  scale  (3,437.74677  nmi  x 
6076  ft/nmi  X  12  in./ft  x  2,000,000).  This  value  translates  to  3,183.293006  mm: 

RAD2*^  =  Earth  radius  @  1:2M 
=  3,183.293006  mm 

CIR^^  =  Earth  circumference  @  1:2M 
=  2  7t  X  RAD^^* 

=  2  X  3.141592654  x  3,183.293006  mm 
=  20,001.21985  mm 

Since  each  row  of  TS  segments  completely  surrounds  the  spherical  earth,  the  number  of  rows 
of  segments  is  calculated  from  the  distance  along  a  semicircle  of  constant  longitude  running  from 
one  pole  to  the  other: 

#  rows  =  (CIR^^/2)/height  of  1  segment  of  paper  chart 
=  (20,001.21985  mm2)/5 1.04373385  mm 
=  195.9223822 

The  number  of  rows  is  truncated  to  195  to  avoid  the  difficulties  in  trying  to  display  a  fractional 
number  of  rows.  Conversion  from  the  data  spacing  in  the  segments  to  the  data  spacing  for  the  equal 
area  display  is  obtained  by  resampling  the  segment  pixel  data  to  the  spacing  of  the  display  pixels. 
This  vertical  resampling  is  called  the  latitudinal  screen  correction  factor  (LSC): 

LSC  =  195/195.9223822 
=  0.9952921 

For  a  particular  map  scale  (known  as  a  “scale  model”),  and  for  a  given  zone,  the  width  and 
height  in  degrees  of  a  segment  is  constant.  The  number  of  segments  per  row  or  column  of  a  given 
zone  is  proportional  to  the  scale.  For  example,  the  1:1M  scale  model  has  twice  as  many  segmf'nts 
around  the  equator  as  the  1:2M  scale  model. 

Since  the  TS  segments  in  the  1:2M  scale  model  contains  195  rows  and  since  the  latitudinal 
height  of  TS  segments  in  all  zones  is  constant  at  any  given  scale,  the  height  of  a  1 :2M  TS  segment 
is  equal  to  1807195,  or  0.92307692°  of  latitude.  The  1:2M  scale  model  contains  304  columns  in 
the  temperate  zones,  and  each  column  has  constant  width  in  degrees  of  lo  gitude.  Thus,  the  width 
of  a  1:2M  temperate  zone  segment  is  360°/304,  or  1.18421053°  of  longitude.  Likewise,  the  1:2M 
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scale  model  contains  380  columns  of  equal  size  in  the  equatorial  zone.  Thus,  the  width  of  a 
1:2M  equatorial  zone  segment  is  3607380,  or  0.94736842°  of  longitude.  Polar  zone  segments  are 
equal  in  size  to  equatorial  zone  segments. 

All  other  scale  models  differ  from  the  1:2M  scale  model  in  that  the  segment  heights  and  widths 
are  proportional  to  the  ratio  between  the  1:2M  scale’s  measurements  and  those  of  the  scale  c  ' 
interest.  For  example,  the  segment  height  and  widths  for  a  1;1M  scale  model  are  determined  as 
follows: 

height  (1:1M)=((1/2,000,000)/(1/I.000,000))x height  (1;2M) 

=  1/2  X  height  (1:2M) 

=  1/2  X  1807195 
=  0.46153846°  latitude 

equatorial  and  polar  width  (1:1M)  =  ((1/2,000,000)/(1/1,000, 000))  x  equatorial  width  (1:2M) 

=  1/2  X  equatorial  width  (1:2M) 

=  1/2  X  3607380 
=  0.47368421°  longitude 

temperate  width  (1:IM)  =  ((l/2,000,000)/(l/l, 000,000))  x  temperate  width  (1:2M) 

=  1/2  X  temperate  width  (1:2M) 

=  1/2  X  3607304 
=  0.59210526°  longitude 


Pixels 

Each  segment  is  subdivided  into  256  rows  and  256  columns.  In  the  equatorial  and  temperate 
zones,  each  of  these  subdivisions  is  bounded  by  lines  of  constant  latitude  and  longitude,  and  the 
latitudinal  or  longitudinal  length  of  a  subdivision  is  equal  to  1/256  times  the  segment’s  latitudinal 
or  longitudinal  length,  respectively.  This  subdivision  is  the  smallest  entity  in  the  TS  model  and  is 
represented  graphically  as  a  pixel.  More  precisely,  the  pixel  is  the  lower-left  comer  point  of  the 
subdivision,  as  shown  in  Fig.  7.  Each  pixel  of  CAC  image  data  is  assigned  an  RGB  color  value. 

Given  the  segment  sizes  in  Table  6,  and  the  fact  that  a  pixel  is  1/256  of  a  segment,  then  the 
size  of  a  1:2M  pixel  can  be  calculsted  as  follows: 

1:2M  pixel  height  (all  zones)  =  0.92307692°/256  =  3.6  x  10”^  degrees  latitude 

1:2M  pixel  width  (polar  &  equatorial  zones) 

=  0.94736842°/256  =  3.7  x  10“^  degrees  longitude 

1 :2M  pixel  width  (temperate  zones) 

=  1.18421053°/256  =  4.6  X  10~^  degrees  longitude 


The  A  .  :  Tessellated  Spheroid  Map  Projection  System:  A  Comprehensive  Definition 


11 
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Fig.  7  —  Pixel  location  within  TS  segment 


Pixel  sizes  for  other  scales  are  directly  proportional  to  the  1:2M  pixel  sizes  (e.g.,  the  longitudinal 
coverage  of  an  equatorial  pixel  at  the  1:1M  scale  is  exactly  one-half  the  longitudinal  coverage  of 
an  equatorial  pixel  at  the  1:2M  scale). 

The  equivalent  6  and  X  values  for  any  given  pixel  can  be  detennined  as  follows: 

Equatorial  and  (rotated)  polar  zone  pixels  (1:2M  scale):  8  =  6*  +  (row**  x  3.7“  x  10~^) 

X  =  X*  +  (column**  x  3.6“  x  ICT^) 

Temperate  zone  pixels  (1:2M  scale):  6  =  8‘  +  (row**  x  4.6“  x  10"^) 

X  =  X*  +  (column**  x  3.6“  x  10“^) 

The  coordinates  of  the  south-west  comer  of  the  segment  in  which  the  pixel  resides  are  8*  and 
X*,  and  row**  and  column**  are  the  row  and  column  of  the  pixel  within  the  segment,  where  (row**, 
column**)  of  the  lower-left  pixel  in  a  segment  is  (0, 0),  and  (row**,  column**)  of  the  upper-right  pixel 
in  the  segment  is  (255,  255). 


Number  of  Rows  of  TS  Segments 

The  number  of  rows  of  1:2M  segments  in  each  zone,  excluding  overlap,  are  as  follows 
(to  include  overlap,  add  four  rows  to  each  figure:  two  rows  for  the  southern  boundary  overlap,  and 
two  rows  for  the  northern  boundary  overlap): 

#  rows  of  NP  segments  =  (upper  NP  boundary  -  lower  NP  boundary)  x  2/segment  height 

=  (90“  -51.69230769“)  x  2/(180“/195) 

=  83 
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#  rows  of  NT  segments  =  (upper  NT  boundary  -  lower  NT  boundary)/segment  height 

=  (51.69230769°  -  31.38461538°)/(1807195) 

=  22 


#  rows  of  EQ  segments  =  (upper  EQ  boundary  -  lower  EQ  boundary )/segment  height 

=  (3 1 .3 846 1 538°  -  (-  3 1 .3846 1 538°))/(  1 80°/ 1 95) 

=  68 


#  rows  of  ST  segments  =  (upper  ST  boundary  -  lower  ST  boundary )/segment  height 

=  (-  31.38461538°  -(-  51.69230769°))/(180°/195) 

=  22 


#  rows  of  SP  segments  =  (upper  SP  boundary  -  lower  SP  boundary)  x  2/segment  height 

=  (-  51.69230769°  ~  (-  90°))  x  2/(180°/195) 

=  83 


Rows  of  polar  segments  run  parallel  to  the  90°  E-90°  W  meridian.  Since  a  rectangular  grid 
cannot  be  evenly  draped  over  a  spherical  surface,  the  number  of  polar  segment  rows  is  the  number 
of  EQ-sized  segments  that  run  along  the  0°-180°  meridian  in  each  polar  zone.  Near  the  “edges”  of 
the  grid,  far  fewer  segments  will  run  in  this  direction,  as  Fig.  5  illustrates.  The  greatest  number  of 
segments  exists  in  the  row  that  crosses  directly  over  the  pole,  along  the  90°  E-90°  W  meridian. 
Each  polar  zone  overlaps  its  neighboring  temperate  zone  by  four  segments  (two  at  each  zone 
boundary).  Note,  however,  that  the  nature  of  the  polar  zones  requires  that  the  polar  overlap  area  is 
not  always  along  a  row  of  polar  segments.  As  an  observer  follows  along  the  edge  of  a  polar  zone 
overlap  region  (for  example,  in  Fig.  1),  the  point  of  overlap  changes  from  being  along  a  row  (at 
0°  longitude)  to  being  along  a  colunm  (at  90°  E  longitude),  then  back  along  a  row  (at  180°  longitude) 
and  again  along  a  column  (at  90°  W  longitude).  Each  nonpolar  zone  overlap  region  is  always 
parallel  to  a  constant  row  of  segments. 


Number  of  Columns  of  TS  Segments 

The  number  of  colunms  of  TS  segments  in  any  nonpolar  zone  utilizes  the  circular  radius 
defined  from  the  axis  of  rotation  of  the  earth  at  a  so-called  “central  latitude”  (CENTXATj). 
CENTXATj,  is  zone-dependent,  as  indicated  by  the  subscript  Z,  but  it  is  independent  of 
scale.  CENTXATj  is  14.123558°  for  the  EQ  zone  and  ±39.120799°  for  the  NT  and  ST  zones.  At 
CENT_LATj,  the  longitudinal  extent  of  a  segment  corresponds  exactly  with  a  51.04373385  mm 
(2.01  in.)  wide  portion  of  chart  and  the  number  of  columns  is  calculated  as  follows: 

#  columns^  =  (CIR^M  x  cos  (CENTXATj))/5 1.04373385 
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Thus,  the  number  of  EQ,  NT,  and  ST  columns  (at  the  1:2M  scale)  may  be  derived  as  follows; 

#  columns  of  1:2M  EQ  segments 

=  NINT  [(20,001.21985  x  cos  (14. 123558‘’))/5 1.04373385]  =  380 

#  columns  of  1;2M  NT,  ST  segments 

=  NINT  [(20,001.21985  x  cos  (39. 120799“))/5 1.04373385]  =  304, 


where  NINT  rounds  to  the  nearest  integer. 

Like  the  latitudinal  correction,  a  longitudinal  screen  factor  correction  is  required  to  convert 
from  data  spacing  in  segments  to  the  equal  area  display.  The  segments  are  resampled  to  the  spacing 
of  the  display  at  the  latitude  of  the  display,  and  the  resampling  spacing  is  called  the  screen  factor 
(SFj).  At  a  given  latitude  LAT  in  a  given  zone  Z,  the  screen  factor  is  calculated  as  follows; 

SF^  =  cos  (CENT_LAT2)/cos  (LAT). 

Due  to  the  orientation  of  the  polar  zone,  the  number  of  columns  of  polar  segments  is  calculated 
in  the  same  way  as  the  number  of  rows  of  polar  segments,  except  that  the  segment  width  is  used 
instead  of  segment  height: 

#  columns  of  NP  segments  =  (upper  NP  boundary  -  lower  NP  boundary)  x  2ysegment  width 

=  (90“  -  51.69230769“)  x  2/(3607380) 

=  80.87 

#  columns  of  SP  segments  =  (upper  SP  boundary  -  lower  SP  boundary)  x  2/segment  width 

=  (-51 .69230769“  -  (-  90“))  x  2/(360“/380) 

=  80.87 

Columns  of  polar  segments  run  parallel  to  the  0“-180“  meridian.  Again,  since  a  rectangular 
grid  cannot  be  evenly  draped  over  a  spherical  surface,  the  number  of  polar  segment  columns  is  the 
number  of  EQ-sized  segments  that  run  along  the  90“  E-90“  W  meridian  in  each  polar  zone.  Near 
the  “edges”  of  the  grid,  far  fewer  segments  will  run  in  this  direction  (Fig.  5).  The  greatest  number 
of  segments  exists  in  the  column  that  crosses  over  the  pole. 

Calculating  TS  Row  and  Column  from  Latitude  and  Longitude 

Table  7  provides  a  partial  list  of  the  row  boundaries  (in  degrees  of  latitude)  for  each  nonpolar 
zone  at  a  scale  of  1:2M,  and  Table  8  provides  a  partial  list  of  the  column  boundaries 
(in  degrees  of  longitude)  for  these  zones  at  the  same  scale.  Table  9  provides  a  partial  list  of  the 
segment  comer  points  (in  both  polar  and  “rotated  equatorial”  latitude  and  longitude  coordinates)  for 
one-quarter  of  the  north  polar  zone  at  a  scale  of  1:2M.  The  complete  versions  of  these  lists  may 
be  found  in  the  appendix  of  the  CAC  Military  Specification  (Lohrenz  1991). 
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Table  7  —  Partial  List  of  Row  Boundaries  (Including 
Overlaps)  for  Nonpolar  Zones  at  1:2M  Scale.  Note  that 
Each  Upper  Row  Boundary  is  1  Pixel  Less  than  the  Lower 
Boundary  of  the  Next  Row  to  the  North. 


Boundaries  (°lat) 
Lower  Upper 


NT  Zone  (Rows  57  to  32): 

57 

52.61538462 

53.53485577 

56 

51.69230769 

52.61177885 

55 

50.76923077 

51.68870192 

54 

49.84615385 

50.76562500 

• 

• 

• 

• 

• 

• 

35 

32.30769230 

33.22716346 

34 

31.38461538 

32.30408653 

33 

30.46153846 

31.38100961 

32 

29.53846154 

30.45793269 

EQ  Zone  (Rows  35  to  -36) 

35 

32.30769230 

33.22716346 

34 

31.38461538 

32.30408653 

33 

30.46153846 

31.38100961 

32 

29.53846154 

30.45793269 

• 

• 

• 

• 

• 

• 

0 

0.00000000 

0.91947115 

• 

• 

• 

m 

• 

• 

-33 

-30.46153846 

-29.53485577 

-34 

-31.38461538 

-30.45793269 

-35 

-32.30769230 

-31.38100961 

-36 

-33.23076923 

-32.30408653 

ST  Zone  (Rows  -33  to  -58): 

-33 

-30.46153846 

-29.53485577 

-34 

-31.38461538 

-30.45793269 

-35 

-32.30769230 

-31.38100961 

-36 

-33.23076923 

-32.30408653 

• 

• 

• 

• 

• 

• 

-55 

-50.76923077 

-49.84254808 

-56 

-51.69230769 

-50.76562500 

-57 

-52.61538462 

-51.68870192 

-58 

-53.53846154 

-52.61177885 

— Zone  Overlap 


Zone  Overlap 


Zone  Overlap 


Zone  Overlap 


Zone  Overlap 


2k>ne  Overlap 
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Table  8  —  Partial  List  of  Column  Boundaries  for  Nonpolar  Zones  at  1;2M  Scale.  Note  that  Each 
Eastern  Column  Boundary  is  1  Pixel  Less  than  the  Western  Boundary  of  the  Next  Colunm  to  the  East. 


Boundaries  (°lon) 

Column 

West 

East 

EQ  Zone  (Columns 

-190  to  189): 

-190 

-180.00000000 

-179.04893091 

-189 

-179.05263157 

-178.10156249 

-188 

-178.10526315 

-177.15419407 

• 

• 

• 

• 

• 

• 

0 

0.00000000 

0.94366776 

• 

• 

• 

• 

• 

• 

187 

177.15789473 

178.10156249 

188 

178.10526315 

179.04893091 

189 

179.05263157 

179.99629934 

Boundaries  (°lon) 

Column 

West 

East 

NT  and  ST  Zones  (Columns  -152  to  151): 

-152 

-180.00000000 

-178.81116361 

-151 

-178.81578943 

-177.62695308 

-150 

-177.63157890 

-176.44274255 

• 

• 

• 

• 

• 

• 

0 

0.00000000 

1.17958471 

• 

• 

• 

• 

• 

• 

149 

176.44736837 

177.62695308 

150 

177.63157890 

178.81116361 

151 

178.81578943 

179.99537418 

The  following  equations  calculate  the  TS  row  and  column  numbers  for  a  given  nonpolar  latitude 
and  longitude  coordinate  (variables  and  functions  are  defined  following  the  equations): 

f_ROW  =  (LAT  X  FACTOR)/(180.0/NTESS’^) 

ROW  =  TRUNC  (f_ROW) 

f.COL  =  (LON  X  FACTOR)/(360.0/NTESS^) 

COL  =  TRUNC  (f.COL) 


To  handle  row  or  column  rounding  and  truncation  errors,  make  the  following  checks  and 
corrections: 

Negative  latitude  case:  If  ROW  <  0  and  I  ROW  -  fJROW  I  ^  PIX  jn'/2,  then  ROW  =  ROW  -  1 
Positive  latitude  case:  If  ROW  >  0  and  I  ROW  -  fJROW  I  <  PIXJHT/2,  then  ROW  =  ROW  +  1 
Negative  longitude  case:  If  COL  <  0  and  I  COL  -  f.COL  I  PIX.WD/2.  then  COL  =  COL  -  1 


Positive  longitude  case:  If  COL  >  0  and  t  COL  -  f.COL  I  <  PIX.WD/2,  then  COL  =  COL  +  1 


Special  cases:  If  COL  ^  NTESS^  x  FACTOR/2,  then  COL  =  COL  -  (NTESS*^  x  FACTOR) 

If  COL  <  -  NTESS^  X  FACTOR/2,  then  COL  =  COL  +  (NTESS^^  x  FACTOR) 
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Table  9  —  Partial  List  of  Segment  Comer  Points  (in  Both  Polar  and  “Rotated 
Equatorial”  Latitude  and  Longitude  Coordinates)  for  One-Fourth  of  the  North 
Polar  (NP)  Zone  at  1:2M  Scale 


ROW  0 

EQLAT 

EQ  LON 

NP  LAT 

NP  LON 

0.00000000 

180.00000000 

90.00000000 

90.00000000 

0.00000000 

179.05263158 

89.05263158 

90.00000000 

0.00000000 

178.10526316 

88.10526316 

90.00000000 

• 

(Total  of  43  Segments  in  Row  0) 

• 

0.00000000 

142.10526316 

52.10526316 

90.00000000 

0.00000000 

141.15789474 

51.15789474 

90.00000000 

0.00000000 

140.21052632 

50.21052632 

90.00000000 

ROW  I 

EQ  LAT 

EQ  LON 

NP  LAT 

NP  LON 

0.92307692 

180.00000000 

89.07697898 

0.00000000 

0.92307692 

179.05263158 

88.67735224 

45.74027634 

0.92307692 

178.10526316 

87.89246947 

64.01953888 

• 

(Total  of  43  Segments  in  Row  1) 

• 

0.92307692 

142.10526316 

52.09571213 

88.49735260 

0.92307692 

141.15789474 

51.14866271 

88.52839661 

0.92307692 

140.21052632 

50.20160046 

88.55780029 

• 

• 

• 

• 

• 

• 

• 

# 

t 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

ROW  42 

EQLAT 

EQ  LON 

NP  LAT 

NP  LON 

38.76923077 

180.00000000 

51.23076897 

0.00000000 

38.76923077 

179.05263158 

51.22101818 

1.17936492 

38.76923077 

178.10526316 

51.19178086 

2.35740948 

• 

(Total  of  13  Segments  in  Row  42) 

• 

3S.16923(m 

170.52631579 

50.26781440 

11.58180809 

38.76923077 

169.57894737 

50.06861610 

12.69224834 

38.76923077 

168.63157895 

49.85160910 

13.78975296 

ROW  43 

EQLAT 

EQ  LON 

NPLAT 

NP  LON 

39.69230769 

180.00000000 

50.30769125 

0.00000000 

39.69230769 

179.05263158 

50.29825586 

1.14121962 

39.69230769 

178.10526316 

50.26996350 

2.28122282 

39.69230769 

177.15789474 

50.22285546 

3.41880035 

39.69230769 

176.21052632 

50.15700025 

4.55275583 

39.69230769 

175.26315790 

50.07249310 

5.68191433 

39.69230769 

174.31578947 

49.96945530 

6.80512571 

39.69230769 

173.36842105 

49.84803339 

7.92127323 
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Variables  are  defined  as  follows: 

ROW,  COL  Integer  row  and  column  number  of  the  TS  segment 

f_ROW,  f_COL  Floating  point  versions  of  ROW  and  COL 

LAT,  LON  Latitude  and  Longitude  of  the  SW  comer  of  the  TS  segment 

FACTOR  Scaling  factor,  which  is  calculated  as  (2,000, 000/(inverse  of  scale)): 

FACTOR  =  1  for  1:2M  scale 
FACTOR  =  2  for  1:1M  scale 
FACTOR  =4  for  1:500k  scale 
FACTOR  =  8  for  1:250k  scale 
FACTOR  =  20  for  1:100k  scale 
FACTOR  =  40  for  1:50k  scale 

NTESS**  Number  of  row-based  tessellations  per  180°  of  latitude  in  the  1;2M 

scale  model: 

NTESS*  =  195  for  all  zones 

NTESS^  Number  of  colunm-based  tessellations  per  360°  of  longitude  in  the  1 :2M 

scale  model: 

NTESS^  =  380  for  Equatorial  and  Polar  zones 
NTESS^  =  304  for  Temperate  zones 

PIX_HT,  PIX_WD  Height  and  width  of  a  pixel  at  a  given  zone  and  scale  (1/256  of  the 

segment  height  or  width,  respectively,  as  defined  in  an  earlier  section 
of  this  report) 

TRUNC  Function  to  truncate  floating  point  number  to  integer;  e.g.: 

TRUNC  (38.7)  =  38,  and  TRUNC  (-26.5)  =  -26 

Ixl  Absolute  value  of  x;  e.g.,  121  =  2,  and  1-21  =  2 

To  calculate  row  and  column  numbers  from  a  polar  latitude  and  longitude,  first  rotate  the  polar 
coordinates  to  equatorial  coordinates  (the  rotation  equations  are  given  in  the  “TS  Polar  Rotations” 
section).  Then,  convert  the  rotated  equatorial  coordinates  (LAT,  LON)  to  TS  ROW,  COL  as  shown 
above  for  nonpolar  data. 
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Calculating  Latitude  and  Longitude  Coordinates  for  TS  Segment 
at  Given  Row  and  Column  Location 

The  following  equations  calculate  the  latitude  and  longitude  of  the  southwest  comer  of  a 
nonpolar  TS  segment  that  is  located  at  a  given  row  and  column: 

LAT  =  (float  (ROW)  x  (180.0/NTESS’^))/FACTOR 

LON  =  (COL  X  (360.0/NTESS‘^))/FACTOR 

Variables  are  defined  in  the  previous  section.  To  calculate  the  polar  latitude  and  longitude  of 
a  TS  row  and  column  number,  first  use  the  above  equations  to  calculate  the  rotated  equatorial 
latitude  and  longitude  (LAT  and  LON)  from  ROW  and  COL.  Then  rotate  LAT  and  LON  to  the 
polar  coordinates  using  the  rotation  equations  given  in  the  following  section  of  this  report. 


TS  POLAR  ROTATIONS 

Rotation  from  Polar  Coordinates  to  Equatorial  Coordinates 

Given  a  pair  of  polar  zone  coordinates  (0^,  X^),  the  rotated  coordinates  (0^,  X^)  of  the 
corresponding  point  in  the  equatorial  zone  are  as  follows: 


arcsin  (cos(0^)  x  cos(X^)]  for  0^  ^  90" 
0"  for  o'*  =  90" 


North  Poiar.  X*  = 


South  Polar  X*  = 


1  arctan 

-cosjo^)  X  sin(X^) 

sin(0^) 

il80" 

1  arctan 

-cos(0^  X  sin(X^) 

sin(<P) 

(180" 

+  180"  for0^?t9O" 
for  0^=90" 

for  o'* 90" 
for  o'*  =  90" 


Figure  8  illustrates  this  transformation  for  both  poles.  One  more  equation  is  included  for 
completeness,  although  it  is  not  applicable  here  (since  the  polar  zone  only  extends  from  ±90"  to 
±51.7"  in  latitude): 


-90"forX^<0" 
90"for  X^^O" 


Special  Case,  if  0^  =  0"  then  X*  = 


Rotation  from  Equatorial  Coordinates  to  Polar  Coordinates 

Given  a  pair  of  rotated  equatorial  zone  coordinates  (8^,  X^),  the  coordinates  (8^,  X^)  of  the 
#  corresponding  point  in  the  polar  zone  are  as  follows: 

North  Polar:  8^  =  arcsin  [cos  (8*)  x  cos  (X*  - 180“)] 


for8*>0“ 

for8*<0“ 

for  (8*«  0“.X*>180“) 
for(8*  =  0“,X*S180“) 


Lohrenz,  Trmchard,  Myrick,  Wisckow,  and  Riedlinger 


20 


South  Polar:  6^  =  arcsin  [  -cos  (0^) 


f 


aictan 


arctan 


+  IgO" 


X  cos  (X*)J 
for0*>O“ 
fore*<0‘* 


90“ 

-90“ 


for(e*  =  0“,X^S0“) 
for(e*  =  0“,X^<0“) 


DOWNSAMPLING  FROM  ARC  TO  TS 

Downsampling  data  from  the  ARC  model  into  TS  reduces  the  image  resolution  by  approximately 
50%:  from  10  to  approximately  S  pixels  per  linear  millimeter  (precisely  254  to  128  ppi).  Conversely, 
downsampling  increases  the  nominal  pixel  spacing  from  100  pm  to  approximately  200  pm,  in  both 
the  east-west  and  north-south  directions.  Thus,  the  amount  of  storage  that  is  required  to  store  a 
given  section  of  chart  data  is  reduced  by  a  factor  of  approximately  4:1  (2:1  in  the  X  direction,  and 
2: 1  in  the  Y).  After  this  transformation,  one  segment  (256  x  256  pixels)  of  TS  data  consists  of 
196,608  bytes;  prior  to  the  conversion  from  ARC  to  TS,  an  equivalent  section  of  map  data  consisted 
of  more  than  780,000  bytes. 

Map  data  are  downsampled  from  ARC  to  TS  with  a  neighborhood  averaging  function.  If  the 
TS  pixel  grid  is  “overlaid”  onto  the  ARC  grid,  there  would  usually  be  four  neighboring  ARC  pixels 
for  each  TS  pixel.  (An  exception  may  occur  at  the  edge  of  an  area,  in  which  case  there  might  be 
less  than  four  neighboring  ARC  pixels).  As  shown  in  Fig.  9,  typically  four  ARC  pixels  (Al,  A2, 


THREE  EXAMPt£S  OF  OOWNSAMPUNG  4  ARC  PIXELS  (At,  A2,  A3,  M)  irm}  m  PIXEL  (T): 

Al . A2  At .  A2  Al - A2 


R,  6,  AND  B  VALUES  OF  TS  PIXEL  -  AVERAGE  OF  TIC  4  ARC  PIXELS' R.  G.  AND  B  VALUES.  RESPECTimV: 

V<AV«q*A3q  +  A4o)/4 
T,.(A1,4A2,*A3b+A4^/4 


Fig.  9  —  Downsampling  four  ARC  pixels  into  one  TS  pixel 
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A3,  and  A4)  define  a  box  around  the  TS  pixel  (T),  but  T  is  not  necessarily  in  the  exact  center  of 
this  box.  To  assign  an  appropriate  RGB  value  to  T,  the  red.  green,  and  blue  components  of  each 
of  the  neighboring  ARC  pixels  are  averaged,  as  shown  in  the  equations  at  the  bottom  of  Fig.  9. 
A  more  general  equation  for  this  operation  is  as  follows: 

r,=  ^  Z  An,  , 

"  n=l 

where  N  =  number  of  ARC  pixels  (up  to  4)  that  neighbor  the  TS  pixel  and  c  =  color  component 
(R,  G,  or  B) 

No  consideration  is  given  to  the  distance  between  the  TS  pixel  and  each  of  the  ARC  pixels  in 
this  function.  We  investigated  the  effects  of  incorporating  a  distance  factor  into  the  equation,  so 
that  the  RGB  value  of  an  ARC  pixel  that  is  closer  to  the  outplit  TS  pixel  would  be  weighted  more 
heavily  than  the  RGB  value  of  an  ARC  pixel  that  was  farther  away.  However,  no  improvement  was 
observed  in  the  quality  of  the  output  image,  and  the  algorithm’s  execution  speed  was  somewhat 
degraded,  so  the  weight  factor  was  not  included  in  the  final  algorithm.  A  follow-up  report  will 
present  this  downsampling  process  in  greater  detail. 


TS  PROJECTION  DISTORTION 

While  there  is  little,  if  any,  observable  distortion  of  CAC  features  in  the  TS  projection,  slight 
distortion  does  occur  during  the  equirectangular  projection  and  display  of  CAC  data  from  scene 
memory.  The  least  amount  of  distortion  exists  at  the  center  of  a  given  display  of  data  (in  the  digital 
moving  map  systems,  that  is  at  the  center  of  a  127  mm^  (5  in.^)  section  of  displayed  chart  data), 
and  distortion  increases  in  the  east-west  direction  as  the  edges  of  the  current  display  are  approached. 
The  variation  in  pixel  spacing  causes  the  distortion  to  accumulate  from  the  meridian  that  passes 

through  the  center  of  the  display.  The 
amount  of  distortion  is  greatest  in  the 
smallest  scales  (i.e.,  worst  in  the  1:2M 
scale)  due  to  the  greater  geographic 
coverage  per  display.  Figure  10  illustrates 
the  distribution  of  projection  distortion 
for  a  127  mm  x  127  mm  display  in  a 
north-up  orientation.  The  greatest  distor¬ 
tion  at  any  given  scale  occurs  at  each  of 
the  four  comers  of  the  display. 

The  maximum  possible  distortion  in 
the  TS  model  occurs  at  the  points  at  which 
the  latitude  is  farthest  from  the  equator 
(i.e.,  at  the  edge  of  the  temperate  zones’ 
pole-most  overlap  regions)  and,  as  pre¬ 
viously  mentioned,  at  the  largest  display 
coverage  (i.e.,  the  smallest  scale,  which 
is  1:2M).  Table  10  gives  the  maximum 
Fig.  10  —  Distribution  of  projection  distortion  distortion  for  each  nonpolar  TS  zone  at 


DISPLAY  DISTORTION  IS  CAUSED  BY 
USING  A  SINaE  lATITUOtNAL  CORREC- 
TKDN  FOR  THE  OITIRE  DISPLAY. 


w 

n  PIXEL  OETORTION 

/ 


484  X  484  0ISPLAY 


ZERO  DISTORTION 
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Table  10  —  Worst-Case  Distortions  (in  Overlap  Segments)  for  Nonpolar 

TS  2k>nes 


Display 

Coverage 

Map 

Latitude 

Distortions 

Distortions 

Scale 

(degrees) 

(pixels) 

(meters) 

(nmi) 

Equatorial  Zone 

1;50.(X)0 

±31.431 

±5.5691  X  10"2 

±5.5256  X  10-' 

±2.9836  X  10-^ 

1:100,000 

±31.477 

±1.1149x  10~‘ 

±2.2124 

±1.1946  X  10-3 

1:250,000 

±31.615 

±2.7956  X  10-‘ 

±1.3869  X  lO' 

±7.4885  X  10-3 

1:500,000 

±31.846 

±5.6186  X  10'* 

±5.5747  X  lO' 

±3.0101  X  10-2 

1:1,000,000 

±32.308 

±1.1347 

±2.2516  X  10^ 

±1.2158  X  10-' 

1:2,000,000 

±33.231 

±2.3129 

±9.1795  X  10^ 

±4.9565  X  10-' 

Temperate  Zones 

1:50,000 

±51.738 

±1.0484  X  10-' 

±1.0402 

±5.6166  X  10*^ 

1:100,000 

±51.785 

±2.0978  X  10"' 

±4.1628 

±2.2477  X  10-3 

1:250,000 

±51.923 

±5.2519  X  10-* 

±2.6054  X  10' 

±1.4068  X  10-2 

1:500,000 

±52.154 

±1.0529 

±1.0446  X  10^ 

±5.6406  X  10-2 

1:1,000,000 

±52.615 

±2.1156 

±4.1982  X  10^ 

±2.2668  X  10-' 

1:2,000,000 

±53.538 

±4.2703 

±1.6948  X  10^ 

±9.1510x  10-' 

Table  11  —  Worst-Case  Distortions  (in  Overlap  Segments) 
for  Polar  TS  Zones 


Rotated 

Polar 

Equatorial 

Display 

Coverage 

Map 

Latitude 

Latitude 

Distortions 

Distortions 

Scale 

(degrees) 

(degrees) 

(pixels) 

(meters) 

(nmi) 

1:50,000 

±51.646 

±38.354 

±6.6275  X  10-2 

±6.5757  X  10-' 

±3.5506  X  10-* 

1:100,000 

±51.600 

±38.400 

±1.3265  X  10"' 

±2.6323 

±1.4213  X  10-3 

1:250,000 

±51.462 

±38.538 

±3.3238  X  10-' 

±1.6489  X  lO' 

±8.9036  X  10-3 

1:500.000 

±51.231 

±38.769 

±6.6729  X  10-' 

±6.6208  X  lO' 

±3.5749  X  10-2 

1:1,000,000 

±50.769 

±39.231 

±1.3446 

±2.6683  X  102 

±1.4407  X  10-' 

1:2,000,000 

±49.846 

±40.154 

±2.7292 

±1.0831  X  103 

±5.8485  X  10-' 

each  display  scale;  Table  11  gives  the  maximum  distortion  for  the  polar  zones  at  each  scale.  The 
derivation  for  these  distortion  values  is  illustrated  in  Figs.  1 1  and  12,  and  the  necessary  variables 
and  constants  are  defined  in  Table  12.  Figure  11a  provides  a  partial  model  of  a  sphere,  and  the 
equations  that  yield  the  precise  width  (before  display)  of  a  segment  at  a  given  latitude  in  the 
equatorial  zone  at  a  scale  of  1:2M.  Figure  11b  provides  a  similar  model  for  the  polar  zone.  Polar 
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R 

=  3,183293  mm  (01:2M  scale) 

R* 

=  Rxsin(90-8) 

=  3,1832%  mm  X  sin  (90-4) 

CIR* 

=2xfW 

= 2x  [3,1 83293  mm  X  sin  (90-4)] 

=20,001 22  mm  X  sin  (90-4) 

W_SEG 

=  CIR^/380  columns  (EQzone) 

=  20,00122  mmx  sin  (90-4)/380 
=  52.635  mm  X  sin  (90-4) 

OUTER  EDGE  Of  POLAR  Z0»€ 
AFTER  ROTATION  OVER 
EQUATORIAL  REGION 


ROTATED  POLAR  ZONE  POINTS  WTO 
GREATEST  01ST0flTX)N  VALUES 
|X  =  0*  FOR  SOUTH  POLAR  TOTATTON. 
180*  FOR  NORTH  POLAR  ROTATION) 


Fig.  11  —  TS  projection  distortion  models:  (a)  nonpolar  zones;  (b)  polar  zones. 


33231*  N 


32J0e*N 


31285*  N 


Ing.  12  —  Display  distortion  (exaggerated)  of  four 
TS  segments  (1:2M  scale,  EQ  zone). 


A  .  44.486  mm  -  44.027  mm 

•  42887x10 ''mm 

•  22129  pixala  (0 128  pNcIt/in.  and  25.4  mnVIn.) 

•  9.1785  xIo’mIO  10^nVmmx2xiaPici8t) 

•  4.9565  X 10-’  nmi  (O  10'*  kmAn.  62137  x  lOr’  mMon. 

and  62806  xIO-'nmMni) 
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Table  12  —  Parameters  Required  in  Calculating  Distortion  for  the  TS  Model 

R  s  Radius  of  earth  at  the  equator,  at  1 :2M  scale 

=  3,183.293006  mm 

=  Radius  of  earth  at  latitude  at  1:2M  scale 
=  R  X  sin  (90  - 

CIR^  =  Circumference  of  earth  at  latitude  at  1:2M  scale 
=  2Jt  X  R^ 

=  2  X  3.1416  X  3,183.293006  mm  x  sin  (90-^) 

=  20,001.21985  mm  x  sin  (90-4) 

NCOL  =  Number  of  TS  columns  at  1:2M  scale 
=  380  for  EQ,  NP,  and  SP  zones 
=  304  for  NT  and  ST  zones 

SCALE  =  Inverse  of  the  data  scale  (e.g.,  if  data  scale  is  1:50,000,  then  SCALE  =  50,000). 

SEG_HT  =  Height  of  1  segment  (‘’latitude);  constant  for  all  zones  within  a  scale,  but  scale-dependent. 

MIN_LAT  =  For  polar  zones  only;  Polar  zone  latitude  that  is  farthest  from  the  pole. 

MAX_LAT  =  Nonpolar  zones:  Latitude  that  is  farthest  from  the  equator  in  the  given  zone. 

Polar  zones:  Rotated  equatorial  equivalent  of  polar  MIN.LAT  latitude  at  0°  longitude. 

W_SEG1  =  Width  (in  mm,  at  1;2M  scale)  of  the  edge  of  a  segment  at  MAX.LAT.  In  the  model  that 
is  used  for  maximum  distortion  calculations,  this  is  the  top  edge  of  the  display: 

_CirMAX_laT/nc0L 

For  EQ,  NP,  and  SP  zones: 

=  [20,001.21985  mm  x  sin  [(90-  MAX_LAT)]/380 
=  52.63478908  ram  x  sin  (90-  MAX_LAT) 

For  NT  and  ST  zones: 

=  [20,001.21985  mm  x  sin  (90  -  MAX_LAT)]/304 
=  65.79345635  mm  x  sin  (90-  MAX.LAT) 

W_SEG2  =  Width  (in  millimeters,  at  1:2M  scale)  of  the  edge  of  a  segment  that  is  1  segment  lower 
than  MAX.LAT.  In  the  model  that  is  used  for  maximum  distortion  calculations,  this  is 
the  center  of  the  di^lay: 

=  CIR<“'^-*^'^'®^-’^/NCOL 
For  EQ,  NP,  and  SP  zones: 

=  (20,001.21985  mm  x  sin  [90  -  (MAX.LAT-SEG_HT)J)/380 
=  52.63478908  mm  x  sin  [90-  (MAXJLAT-SEG.HT)] 

For  NT  and  ST  zones: 

=  (20,001.21985  nun  X  sin  [90  -  (MAX_LAT-SEG_HT)])/304 
=  65.79348635  mm  x  sin  [90-  (MAXJLAT-SEG.HT)] 

=  Maximum  display  distortion  (at  1:?M  scale)  in  millimeters: 

=1  W.SEGl  -  W.SEG2  I 

A**  =  Maximum  display  distortion  in  pixels  (9  exactly  128  pixels/in.  and  25.4  mm/in.): 

=  A”®  X  5.03937008  pixels/mm 

A®  =  Maximum  coverage  distortion  (converted  to  "real  world”  1:1  scale)  in  meters: 

=  A®®  X  SCALE  x  0.001  m/mm 

A"  =  Maximum  coverage  distortion  in  nautical  miles: 

=  A®  X  (1  km/1000  m)  X  (1  mi/1.609344  km)  x  (0.86897624  nmi/mi) 

=  A®  X  0.000539%  nmi/m 
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zone  distortions  were  calculated  from  the  rotated  equatorial  coordinates  of  the  points  that  are 
farthest  from  the  equator  on  the  rotated  equatorial  grid  (i.e.,  at  the  edge  of  the  overlap  region  and 
lying  on  the  0°-180®  meridian).  Figure  12  illustrates  the  distortion  of  four  TS  segments:  the  dashed 
lines  indicate  how  the  segments  would  appear  if  there  were  no  distortion,  and  the  solid  lines 
indicate  how  the  segments  are  actually  displayed.  (Note  that  this  figure  significantly  exaggerates 
the  distortion  in  order  to  illustrate  its  effect).  The  maximum  display  distortion  A  is  calculated  as  the 
difference  between  the  ideal  segment  width  (which  is  the  width  of  the  segment  at  the  center  of 
the  display)  and  the  predistorted  width  of  the  segment  at  the  edge  of  the  display.  Figure  12  presents  the 
distortion  (in  display  millimeters,  pixels,  and  nautical  miles  on  the  Earth)  for  a  segment  at  the  edge 
of  the  equatorial  zone  at  a  scale  of  1:2M.  The  illustrated  distortion  values  in  pixels,  meters,  and 
nautical  miles  correspond  to  the  distortion  values  in  Table  10  for  this  zone  and  scale. 


APPLYING  TS  TO  DATABASES  OTHER  THAN  CAC 

The  TS  projection  system  was  originally  designed  to  store  CAC  data  around  the  globe  in  a 
seamless  manner  for  use  in  aircraft  mission  planning  systems  and  digital  moving  map  systems. 
However,  the  presentation  of  the  CAC  data  in  TS,  along  with  the  efficient  compression  strategy  that 
is  used,  has  piqued  the  interest  of  several  potential  users  of  other  databases.  Two  in  particular  are 
a  proposed  Compressed  Nautical  Chart  (CNC)  database  and  a  compressed  version  of  the  DMA 
Digital  Land  Mass  System  (DLMS)  database.  The  CNC  would  be  identical  to  the  CAC,  except  that 
the  data  would  consist  of  nautical  charts  instead  of  aeronautical  charts.  The  primary  concern  for 
such  a  database  would  be  scale:  the  maximum  scale  required  for  nautical  chart  data  is,  logically, 
considerably  larger  than  the  maximum  scale  that  is  required  for  aeronautical  chart  data.  DMA  plans 
to  provide  the  uncompressed  nautical  chart  data  in  the  form  of  ADRG,  and  NRL  has  proposed  to 
transform  this  nautical  ADRG  data  into  TS  and  compress  it  in  exactly  the  same  way  that  CAC  is 
processed  now. 

The  compression  of  DLMS  would  be  completely  different  from  the  CAC  and  CNC  compressions, 
since  DLMS  is  fundamentally  different  from  ADRG.  DLMS  is  a  union  of  Digital  Terrain  Elevation 
Data  (DTED)  and  Digital  Feature  Analysis  Data  (DFAD).  DTED  consists  of  elevation  posts  that  are 
equally  spaced  in  latitude  and  longitude.  DFAD  is  a  list  of  point,  linear,  and  areal  features  that 
are  stored  as  single  points,  vectors,  and  polygons.  Neither  dataset  is  stored  in  any  specific  map 
projection,  so  the  “transformation”  of  DLMS  to  TS  merely  consists  of  splitting  the  data  into  discrete 
TS  segments,  according  to  their  latitude  and  longitude  positions.  The  scale  of  DLMS  is  determined 
by  the  DTED  post-spacing  and  is  considered  to  be  closest  to  the  l:2S0k  TS  scale  model. 

The  following  sections  will  discuss  several  issues  that  must  be  addressed  before  TS  can  be 
applied  as  a  standard  projection  and  display  system  for  databases  other  than  CAC.  All  of  these 
issues  are  associated  with  scale. 

Distortion  Due  to  Scale  Incompatibility 

The  model  is  currently  configured  to  handle  any  scale  that  can  be  integrally  divided  by 
1 :2M  (the  base  scale  for  the  TS  projection  system).  For  example,  nautical  chart  data  at  an  original 
scale  of  1:7091  would  have  to  be  resampled  into  a  1:8k  scale  model  before  it  could  be  considered 
compatible  with  the  TS  model  (8,000  is  the  closest  number  to  7091  that  can  be  integrally  divided 
into  2,000,000). 
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Such  resampling  is  not  a  major  problem  for  chart  data  that  will  be  converted  to  a  scale  relatively 
close  to  the  original  chart  scale,  but  some  distortion  problems  could  result  from  resampling  data  to 
a  scale  that  is  not  reasonably  similar  to  the  original.  For  example,  when  data  are  resampled  to  a 
smaller  scale  (e.g.,  from  1:7091  to  1:8000),  each  pixel  in  the  dataset  is  forced  to  represent  a  larger 
area  (in  this  case,  each  resampled  pixel  must  represent  an  area  that  is  approximately  1.13  times 
larger  than  the  area  that  was  represented  by  the  original  pixel).  Likewise,  when  data  are  resampled 
to  a  larger  scale  (for  example,  from  1:8652  to  1:8000),  each  pixel  in  the  resampled  dataset  will 
represent  a  smaller  area  than  a  pixel  in  the  original  dataset  (in  the  given  example,  each  resampled 
pixel  would  represent  an  area  that  is  approximately  0.92  times  the  size  of  the  original  pixel).  In 
either  case,  some  distortion  can  be  expected  in  the  geographic  integrity  of  the  resampled  dataset, 
and  the  amount  of  distortion  would  be  directly  proportional  to  the  degree  of  resampling  that  must 
be  done. 

TS  Segment  Row  Limits 

TS  segment  filenames  are  derived  from  the  segment’s  row  and  column  values  and  the  TS  zone 
in  which  it  resides.  As  described  here,  the  maximum  number  of  rows  (covering  the  earth  ellipsoid) 
supported  by  the  TS  model  is  9000,  and  there  are  19S  rows  of  TS  segments  at  the  1:2M  base  scale. 
The  number  of  rows  in  any  other  scale  is  directly  proportional  to  the  ratio  of  the  base  scale  to  the 
scale  in  question.  For  example,  the  1:1M  scale  has  twice  as  many  rows  as  the  1:2M  scale.  Thus, 
the  largest  scale  that  can  be  supported  in  the  current  TS  model  can  be  calculated  as  follows: 

largest_scale  =  2,000,000  x  (195/9000)  =  43,333.33 

Since  the  scale  in  question  must  also  be  an  integral  divisor  of  the  base  scale  (which  1:43.333 
is  not)  the  largest  scale  that  is  currently  supported  by  the  TS  model  is  1:50k. 

However,  CNC  data  will  include  nautical  charts  at  scales  as  large  as  1:8k.  Several  methods  are 
being  investigated  that  would  allow  these  larger  scales  to  be  represented  in  the  TS  model  without 
having  to  significantly  modify  the  specification  for  TS  or  for  the  CAC/CNC  products.  For  example, 
the  segment  filenames  could  be  encoded  in  other  than  base'^  for  scales  larger  than  1 :50k.  Hex 
(base*®)  numbers,  base^  (which  would  use  the  letters  A  through  Z  instead  of  numerals  0  through 
9),  or  base^®  (which  would  use  the  numerals  0  through  9  and  letters  A  through  Z)  are  all  being 
considered  for  scales  larger  than  1:50k.  Base***  would  be  maintained  for  scales  smaller  than  1:50k, 
so  that  existing  systems  that  use  the  CAC  database  (but  which  do  not  plan  to  use  the  CNC  database) 
would  not  be  affected.  Systems  that  do  plan  to  use  CNC  would  have  to  include  a  check  for  scale, 
and  then  decode  the  filenames  according  to  the  appropriate  base  system.  A  future  report  will 
identify  which  base  is  chosen  for  these  larger  scales. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  TS  model  (formerly  referred  to  as  TS  Model  IV,  or  TS-4)  provides  a  seamless,  global 
framework  in  which  to  store  scanned  chart  data  and,  potentially,  other  data  types  as  well.  TS 
sections  the  data  into  latitudinally  based  zones  and  then  subdivides  the  zones  into  tessellations,  or 
segments,  of  256  x  256  pixels.  The  TS  model  was  originally  conceived  for  the  CAC  database, 
which  includes  six  scale  models  of  chart  data  that  range  from  1:50k  to  1:2M.  All  TS  parameters 
are  based  on  the  1:2M  scale  model,  and  any  additional  scales  that  are  to  utilize  TS  must  foe  an 
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integral  divisor  of  2,000,000.  In  addition,  the  TS  model  currently  stipulates  that  segment  filenames 
are  dependent  on  a  maximum  of  9000  TS  segment  rows;  however,  scales  that  are  larger  than 
approximately  1:50k  have  more  than  9000  rows.  These  restrictions  do  not  necessarily  prohibit 
additional  scales  from  being  included  in  the  TS  model,  but  they  must  be  resolved  before  any  new 
scales  are  added.  Several  potential  solutions,  such  as  the  utilization  of  some  base  other  than  base‘s 
to  encode  the  segment  filenames,  were  suggested.  A  future  report  will  identify  the  base  chosen  for 
these  larger  scales.  Another  follow-up  report  will  present  the  downsampling  process  (interpolation 
of  TS  pixels  from  ARC)  in  greater  detail. 
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